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ABSTRACT 

Helioseismology has revealed many details of solar differential rotation and also its time variation, 
known as torsional oscillations. So far there is no generally accepted theoretical explanation for 
torsional oscillations, even though a close relation to the solar activity cycle is evident. On the 
theoretical side non-kinematic dynamo models (including the Lorentz force feedback on differential 
rotation) have been used to explain torsional oscillations. In this paper we use a slightly different 
approach by forcing torsional oscillations in a mean field differential rotation model. Our aim is not a 
fully self-consistent model but rather to point out a few general properties of torsional oscillations and 
their possible origin that are independent from a particular dynamo model. We find that the poleward 
propagating high latitude branch of the torsional oscillations can be explained as a response of the 
coupled differential rotation / meridional flow system to periodic forcing in mid-latitudes, of either 
mechanical (Lorentz force) or thermal nature. The speed of the poleward propagation sets constraints 
on the value of the turbulent viscosity in the solar convection zone to be less than 3 x 10®m^s~^. 
We also show that the equatorward propagating low latitude branch is very unlikely a consequence 
of mechanical forcing (Lorentz force) alone, but rather of thermal origin due to the Taylor-Proudman 
theorem. 

Subject headings: Sun: interior — rotation — helioseismology — dynamo 



1. INTRODUCTION 

Solar torsional oscillations have been known to exist 
for more than two decades. iHoward fc Labont^ l|198Ciri 
presented the first observations of torsional oscillations 
using Mt. Wilson Doppler measurements and pointed 
out the 11 year periodicity and the relation to the so- 

\ lar cycle. These early observations showed only the 
equatorward propagating branch at low latitudes. The 
high latitude branch (above 60°), which is in amplitude 
at least twice as strong as the equatorward propagat- 
ing branch, was foun d more recently through helioseis- 
mic me asurements byl Toomre et-alJ (l2QQfll) ■ LHi)WB_eii-an 

; l[200il. IAntiafc Ra,siin2nnn. Vorontsov et a,].l ^20021). 
and 'Howe et al.l l|2005(l . These inversions also show that 
the high latitude signal penetrates almost all the way 
to the base of the convection zone. The depth penetra- 
tion of the low latitude signal is more uncertain due to 
the lower amplitude, which is comparable to the uncer- 
tainties of the inversion methods in the lower half of the 
convection zone. The most interesting feature of the low 
latitude branch is an inclination of the phase with respect 
to the rotation axis by about 25° . This inclination is very 
close to the inclin ations of the isorot ation contours of Vl 
as pointed out bv lHowe et alJ l|2004D . 

On the theoretical side a variety of explanations have 
been proposed: macroscopic Lorentz force feedback, mi- 
croscopic Lorentz force feedback, and thermal forcing. 

The idea of macroscopic Lorentz force feedback (com- 
puted from the large scale magnetic mean field of the so- 
lar dynamo) wa, s origin ally proposed by|[£chussler (1981) 
and Yoshi mural l|1981|) and has bee n incorporat e d into 
dyna.mo rn odels more recently by iCovas et al.l l|2000l 
120041 l2005|l . While these models address the non-linear 

*The National Center for Atmospheric Research is sponsored by 
the National Science Foundation 
Electronic address: rempel@hao .ucar.edu| 



Lorentz force feedback using a simplified equations of 
motion (c onsidering only the l ongitudin al component), 
models bv I.TenninersI (llflQ,"^ and iRemnell |2006) consider 
the Lorentz for ce feedback al s o in t he meridio n al pla ne. 
The model of Re mpel et all (|2005|l . iRempell 1I2OOI is 
along the lines of the aA-models by 'Brandenburg^^^U 
li^gO, 1991, 1992), Moss ct al. (1995), and Muhli et alJ 
(|1995l|) (coupling mean field models for differential rota- 
tion, meridional fiow and magnetic field evolution), but 
puts more emphasis on the role of the meridio nal flow 
leadin g to a flux-transport dynamo (see e.g. iDikpati 
l|2005(l for a recent review on the development of flux- 
transport dynamos). 

Microscopic Lorentz force feedback (quenching 
of turbulent transport processes driving differen- 
tial rotation 'A- quenching') has been addressed by 
Kitchatinov & Pioin (1998), Kitchatinov ct al. (1999), 
and Kiiker et al . (1999).^^ 

Very recently iSoruit (|2003|) proposed a thermal origin 
of the low latitude branch of torsional oscillations, driven 
through enhanced radiative losses in the active region 
belt. This theory also predicts an inflow i nto the active 
region belt, which h as been observed by Komm et ajj 
CHI, and lZhan Wl^ichev (2004) . 

The inve s tigati on presented here is based on the model 
of iRempeil (|2006ft : however we take a different approach 
by decoupling the excitation of torsional oscillations from 
a detailed dynamo model. We use the differential rota- 
tion model of .Rempel (2QQ5. ) and force torsional oscilla- 
tions through mechanical and thermal perturbations to 
address the question of what type of forcing is required 
to get a response comparable to the observed torsional 
oscillations. 

We emphasize that this is not a fully self-consistent 
explanation of torsional oscillations since we do not in- 
corporate a full dynamo model also considering the non- 
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Fig. 1. — a) Differential rotation (solid contours correspond to faster rotation, dashed contours to slower rotation than the core of the 
sun); b) Stream function of meridional flow (solid contours correspond to a counter clockwise flow). The amplitude of differential rotation 
is 30% of the core rotation rate, the amplitude of the meridional flow is flms"-^ (2ms~-'^) at the top (bottom) of the convection zone. 



linear feedback of zonal and meridional flow variations 
on the magnetic field evolution. We have published such 
a self-consistent non-linear dynamo model (based on a 
flux-transport dynamo and the same differential rotation 
model we use in this paper) in^smocl (2006). We found 
that many features of torsional oscillations discussed in 
that paper are more general, meaning independent from 
a particular dynamo model. The goal of this paper is 
to point out these general properties that are strongly 
related to influence of the Taylor-Proudman theorem on 
amplitude and phase of torsional oscillations. 

2. MODEL 

Our theoretical analysis i s based on the differential ro- 
tation model developed bv iRempell l)2005j) . This model 
is an axisymmetric mean field model incorporating a 
parametrization of turbulent angular momentum trans- 
po rt 'A-effect' (K itcliatinov & Riidiecr 1993). We refer 
to IRempell l)2005|) for details of this model. For this in- 
vestiga tion we use a r eference model close to the model 
used in lRempeil l)2006|) . with the minor difference of using 
a larger value of S = —3 x 10~^ for the superadiabatic- 
ity in the overshoot region (leading to more solar like 
differential rotation in terms of inclination of the isoro- 
tation contours). Figure^shows the differential rotation 
and the meridional flow stream function for the reference 
model. 

Torsional oscillations can be driven in general through 
mechanical forcing (Lorentz force) or thermal forcing 
(pressure imbalances drive geostrophic flows in a rotating 
system). Most investigations so far focused on the role 
of the longitudinal Lorentz force, neglecting azimuthal 
components. This can be partially justified by the fact 
that the toroidal field is significantly stronger than the 
poloidal field in an afl dynamo, leading also to a dom- 
inant component of the Lorentz force in the longitudi- 



nal direction. The azimuthal component could be im- 
portant for explaining (observed) deviations from the 
Taylor-Proudman state through a magnetostrophic bal- 
ance; however, as we show later, this is unlikely to hap- 
pen in the bulk of the convection zone. Thermal forcing is 
very efficient in that reproducing the observed amplitude 
requires only temperature fluctuations of a few tenth of 
a degree, which can be derived from a geostrophic bal- 
ance. Assuming a latitudinal balance between Coriolis 
force and pressure force gives: 



O O 2 . Idp ^lAp ^ p AT 



which leads to 



AT « M {Vlf^rf sin(26l) (^\(- 

\noj \n 



(1) 



(2) 



Here Hq denotes the reference state rotation rate, fii 
the perturbation (torsional oscillation), B is the gas 
constant, fi the mean molecular weight, r the radius, 
and 9 the co-latitude. Adopting solar values of fio — 
2.7 X 10"^ s-\ r = 7 X IQS m, ^ = 0.62, A0 = 10° (width 
of active region belt) and fli/flo ^ 0.005, we get 



AT - 0.25 K 



(3) 



This value is cons istent with the temperature perturba- 
tion of 0.1 to 0.2 K lRemneil (|200<?) imposed in their model 
to drive the low latitude branch of torsional oscillations. 
We emphasize that the total temperature perturbation 
is nearly independent of depth for a fixed torsional oscil- 
lation amplitude, meaning that the relative perturbation 
close to the surface is around 5 x 10~^ (the same value 
was given by Soruit (20(^), but only around 10~^ close 
to the base of the convection zone. 

The following questions will be addressed in this inves- 
tigation: 
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Fig. 2. — Torsional oscillations produced by periodic perturbations at a fixed location. The left panels (a, c, e) show the response to 
a periodic perturbation in the angular velocity (mechanical forcing), the right panels (b, d, f) the response to a periodic perturbation 
in entropy (thermal forcing). Closed contours correspond to positive values (faster rotation), dashed contours to negative values (slower 
rotation). From top to bottom the position of the perturbation in latitude is varied between 60° (a, b), 45° (c, d), and 30° (e, f). In all 
cases the strongest signal is found at the pole, independent from the location of the forcing. The power on the poleward side is always 
higher than on the equatorward side. 



Are torsional oscillations of thermal or mechanical 
origin' 



? 



• Are they a consequence of traveling or periodic per- 
turbations in the solar convection zone? 

We focus our discussion separately on the high lati- 
tude and low latitude branch for reasons that will become 
more evident in the following discussion. 

3. HIGH LATITUDE BRANCH 

The high latitude branch has a larger amplitude than 
the low latitude branch and is propagating poleward, 
while the low latitude branch is clearly following the 
equatorward propagating magnetic activity belt. This 
leads to the fundamental question of whether the pole- 
ward branch is also associated with a poleward propa- 
gating magnetic pattern at the base of the convection 



zone (as it has been seen in many interface dynamos), 
which is not visible at the solar surface. We cannot rule 
out such a possibility here, but we will show that this 
is not conclusive and the poleward propagation can be 
explained as a response of the coupled differential rota- 
tion meridional flow system to a non-traveling periodic 
perturbation in mid-latitudes. To demonstrate this, we 
incorporate in the differential rotation model a mechan- 
ical (in the equation for 17) or thermal (in the entropy 
equation) forcing term of the form 

' .. + Anf{r,e)sm{Lo,t) (4) 



at 

dsi 
'dt 



Asf{r,d) sin(wci) ■ 



(5) 



Here ujc corresponds to a 11 year periodicity (as observed 
for torsional oscillations) , Aq and Ag are the amplitudes 
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Fig. 3. — Influence of viscosity and meridional flow on tlie poleward propagation of zonal flow variations. Panels a) and b) compare to 
models different values of the turbulent viscosity: a) = 3 X 10* ra? and b) ut = 1-2 X 10^ ro? . A higher value of the turbulent 
viscosity leads to a faster poleward propagation of the flow pattern. Both cases show the response to a periodic mechanical forcing at 45° 
latitude. Panel c) shows the same model as panel a), however, the meridional is not allowed to change in response to the zonal flow variation. 
In this case the maximum amplitude coincides with the location of the forcing region and the velocity of the poleward propagation of the 
signal is reduced. Panel d) shows the variation of meridional flow pattern (the reference state meridional flow is subtracted) associated 
with poleward moving zonal flow pattern shown in Panel a). The amplitude of the meridional flow variation is around 0.2ms~^ for a solar 
like zonal flow variation of around 4 nHz at the pole. 



of the forcing, and /(r, 9) is a Gaussian profile of the 
form: 



/(r, 0) = cxp 



Ar 



•exp 



AO 



(6) 



We use for all following models the parameters vq — 
0.85 i?0, Ar = 0.05 i?0, A9 = 0.125, and vary the co- 
latitude between 7r/6, 7r/4, and 7r/3. For the results pre- 
sented here the amplitudes ^ji and are of secondary 
concern (as long as 17i <C flo)- We have chosen values 
leading to torsional oscillations of a few percent of the 
rotation rate. 

Figure 121 shows on the left side models with mechanical 
forcing and on the right side models with thermal forcing. 
From top to bottom, the latitude of the forcing is varied 
between 60°, 45°, and 30°. Independent of the nature 
of the periodic forcing, all models show a strong pole- 
ward propagating pattern and a very weak equatorward 
propagating pattern away from the latitude at which the 
forcing is applied. The equatorward pattern increases 
in amplitude when the location of the forcing is moved 
closer to the equator. The pattern typically shows two 
peaks, one close to the latitude of the forcing and one 
right at the pole (with larger amplitude). Comparing 
the latitudinal extent of the pattern on the polar side of 
the forcing region with observations, a forcing location 
in mid-latitudes seems most reasonable. 



The time for the signal to travel from about 60° to 
the pole is 3 to 4 years, which is close to the observed 
propagation speed. It turns out that the travel time is 
primarily affected by the turbulent viscosity assumed in 
the reference model, which is shown in Figure |21 To this 
end we compare two models with the viscosity values 
of 3 X lO^m^s"^ (panel a) and 1.2 x lO^m^s"^ (panel 
b). We have adjusted in each model the amplitude and 
the direction of the turbulent angular momentum trans- 
port (A-effect) such that the differential rotation and 
meridional flow remain roughly unchanged. The prop- 
agation time for the signal drops from around 3 to 1 
year when the turbulent viscosity is increased by a fac- 
tor of 4. Since solar observations show a time delay of 
3 to 4 years, this sets constraints on the value of the 
turbulent viscosity. In our model we get the best fit for 
values ^ 3 X 10^ m^ s~^. This value is around one order 
of magnitude smaller than mixing-length estimates for 
the solar convection zone, it agrees however with the re- 
sults of Kitchatinov ct al. (1994) taking into account ro- 
tational quenching of turbulent viscosity. Panel c) shows 
the same model as panel a) but the meridional flow is 
not allowed to change in response to the zonal flow vari- 
ation (we solve only the equation for and keep the 
meridional flow fixed). This leads to a significantly dif- 
ferent flow pattern: the maximum amplitude is found 
close to the forcing region and the propagation of the 
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Fig. 4. — Snapshots of the torsional oscillations caused by mechanical forcing (left) and thermal forcing (right). The snapshots correspond 
to t = years in Fig. |^ Torsional oscillations caused by thermal forcing have the tendency to show stronger deviations from the Taylor- 
Proudman state. 



signal toward the pole slows down. In this way the vari- 
ation of the meridional flow has a significant influence on 
the structure of the torsional oscillation pattern. Inter- 
estingly, the meridional flow of the reference state is not 
important, only the induced flow variation. Constructing 
a model with the same value of the turbulent viscosity, 
but a different meridional flow, leads to a nearly identical 
zonal flow pattern (difference in the 1% range). Panel d) 
shows the meridional flow variation associated with the 
poleward moving zonal flow pattern shown in panel a). 
A solar like amplitude of around 4 nHz for the torsional 
oscillation leads to a meridional flow variation of around 
0.2ms-i. 

To our knowledge this is so far the only way to con- 
strain the turbulent viscosity directly by observations, 
without having to relate it to the magnetic diffusivity by 
making assumptions about the magnetic Prandtl num- 
ber. We emphasize that our model assumes a constant 
value throughout the convection zone. In models with 
a depth dependent viscosity the above mentioned value 
should reflect more the values in the lower half of the 
convection zone rather than the surface layers, since they 
contribute more to the angular momentum transport due 
to their larger density. 

From the results shown in Figure [21 it is difflcult to 
judge whether it is possible to distinguish between a me- 
chanical or thermal forcing. Figure ^ shows meridional 
cross-sections of the zonal flow pattern. The most obvi- 
ous difference here is that the mechanical forcing leads 
to torsional oscillations more aligned with the rotation 
axis (Taylor-Proudman theorem), while thermal forcing 
shows a tendency to produce more radially aligned pat- 
terns. Beside aligning structures with the rotation axis, 
the Taylor-Proudman constraint also leads to an almost 
constant amplitude within the convection zone, while in 
the case of thermal forcing the entropy perturbations al- 



low for more variation of the amplitude as function of 
depth. 

In the case of the high latitude branch the observations 
are not detailed enough to distinguish between thermal 
and mechanical forcing (this may be because in high lat- 
itudes the radial direction almost coincides with the axis 
of rotation). We will show in the next section that there 
is a significant difference between mechanical and ther- 
mal forcing effects in low latitudes. 

4. LOW LATITUDE BRANCH 

The most striking feature of the low latitude branch 
is the systematic deviation fro m the Tay lor-Proudman 
state. As first pointed out bv iHowe et al. (2004), the 
lines of constant phase show an inclination of 25° to the 
axis of rotation, similar to the isorotation contours of 
the differential rotation. We present here a few idealized 
experiments to differentiate between possible mechanical 
and thermal forcing. To this end we construct a mechan- 
ical forcing function and a cooling function that has a 25° 
inclination angle with respect to the axis of rotation. For 
reasons of simplicity we apply here a stationary pertur- 
bation for a time interval of a few month to illustrate the 
effect. The forcing functions are given by: 

' ..+An9{r,9)[e-ar)] (7) 



dt 
dsi 

'dt 



+ A,g{r, 9) 



(8) 



with 



g{r, 0) = O.5 



1 -|- tanh 



Ar 



■ exp 



Cir) 



AO 



(9) 



The function C(r) specifies the inclination of the forcing 
with respect to the rotation axis. Assuming that is 
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Fig. 5. — Mechanical forcing of low latitude zonal flow variations. Panels a) and c) show the forced zonal flow considering no change of 
the meridional flow (solving the equation for Oi for the fixed reference state flow), while panels b) and d) show the zonal flow solving the 
full system. Panels a) and b) use a forcing function with = O.75i?0,9o = 80°i c) and d) with rg = 0.95 Rq, 60 = 65°. While panels a) 
and c) show a zonal flow that closely resembles the properties of the forcing function [especially the 25° inclination in panel a)] , the results 
shown in panels b) and d) are strongly influenced by rotation. The flow pattern becomes strongly aligned with the axis of rotation and 
also the amplitude of the zonal flow band with slower rotation is significantly smaller than the faster zonal band. 



the location of the forcing at and A is the inchnation 
angle we get: 



CM = A 



arcsm 



sin(6'o — A) 



(10) 



Here 9q determines the co-latitude of the perturbation 
at the radius rg, where rp marks the depth at which the 
forcing decreases to half of the convection zone value. 
We use in the following discussion common values of 
Ar = 0.05 i?o, A = 25°, and A9 = 0.1 {A0 = 0.05) 
for mechanical (thermal) forcing, respectively. 

To compare results for forcing functions extending 
through the entire convection zone with those from sur- 



face forcing functions, we will use in both cases the 
parameters ro — 0.75 Rq, 60 = 80° (co-latitude); and 
ro = 0.95 i?o, 6*0 = 65°. The amplitudes Aq and As are 
chosen so that applying the forcing for a time interval of 
about 4 months leads to an amplitude of the zonal flow 
of a few uHz. All results shown in the following discus- 
sion show the zonal flow after 4 months of mechanical or 
thermal forcing. 

The forcing function for the mechanical forcing Eq. Q 
is multiplied by the term [0 — C,{r)\ in order to produce 
a perturbation that has a faster zonal flow equatorward 
and slower zonal flow poleward of the region the fore- 
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Fig. 6. — Thermal forcing of low latitude zonal flow variations. The panels a) and c) show the imposed entropy perturbation, the panels b) 
and d) the resulting zonal flows. Panels a) and b) use a forcing function with = 0.75 Rq , Oq = 80° , c) and d) with rg = 0.95 Rq , 6o = 65° . 
Contrary to the results presented in Figure |S| for mechanical forcing, the driven zonal flows show significant deviations from the Taylor- 
Proudman state. Especially in panel b) the imposed 25° angle of the thermal forcing function is clearly visible in the zonal flow field. 
Also both zonal bands (faster and slower one) are visible, with the faster band on the equatorial side being of larger amplitude and more 
confined. 



ing is applied. In the case of the thermal forcing this is 
not required, since the response to a cooling automati- 
cally leads to the formation of two geostrophic flows with 
opposite directions. 

FigureElsummarizes the results obtained with mechan- 
ical forcing, Figure^lthe results from thermal forcing. In 
the case of mechanical forcing the zonal flows do not re- 
semble the properties of the forcing function (especially 
the 25° inclination) unless meridional flow variations are 
completely suppressed as shown in Figure and c) . 
If meridional flow variations are considered, the zonal 
flow pattern shows strong influence of rotation, aligning 



the pattern with the axis of rotation (Taylor-Proudman 
theorem). As a consequence, the resulting flow is en- 
tirely dominated by the influence of rotation and does not 
resemble the properties of the applied forcing function. 
Since the pattern has the tendency to spread parallel to 
the axis of rotation, a perturbation traveling at the base 
of the convection zone (r — 0.72 i?©) from 30° latitude 
toward the equator, would produce close to the surface a 
pattern moving from around 40° to 20° and could there- 
fore not explain the observed low latitude branch of the 
observed torsional oscillations. Even applying the forc- 
ing close to the surface does not lead to a reasonable flow 
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Fig. 7. — Combination of mechanical and thermal forcing; a) 
pure mechanical forcing, b) pure thermal forcing, c) combination 
of both. The mechanical and thermal forcing functions used here 
have a radial alignment and the amplitude of the zonal flow in 
a) and b) is roughly the same. The combination of both forcings 
produces roughly a 25° inclination as observed in the solar case. 



pattern as shown in Figure [SJi). 

This result is not too surprising, since adding a forcing 
in the </> direction does not impact the meridional force 
balance that leads to the Taylor-Proudman state. This 
raises the question if additional magnetic stresses in the 
meridional plane could cause a deviation from the Taylor- 
Proudman state as observed (magnetostrophic balance). 
In a typical afJ-model for the solar dynamo, the poloidal 
field is at least a factor of 100 weaker than the toroidal 
field. Having a toroidal field strength of around 1 T (10 
kG) (Rempcl .2006J, this yields around 0.01 T (100 G) for 
the poloidal field. An estimate similar to Eq. Q leads 
in the case of a magnetostrophic balance to 



(11) 




Here we assumed that the ^-component of the Lorentz 
force can be estimated as Bp / (/io r A9) , where Bp de- 
notes the poloidal field strength. Adopting solar val- 
ues of rio = 2.7 X 10-6 s-i, r = 7 X lO^m, A6I = 10°, 
i^i/ilo ^ 0.005, and we get 



(12) 



where Qhc — 200 kg m^^^ denotes the density at the base 
of the convection zone. These values indicate that a mag- 
netostrophic balance is very unlikely in the bulk of the 
convection zone, but very close to the surface, where the 
density is very low, it cannot be ruled out (for a typical 
solar model we have g « 10~^gibc at r = 0.98 i?© and 
g « lO-'^gbc at r = 0.997 i?©). There are additional 
terms of the Lorentz force in the meridional plane that 
include the much stronger toroidal field: the gradient of 
the magnetic pressure and magnetic tension force arising 
from spherical geometry ~ B'^/r. The latter can lead to 
a prograde iet wit hin the ma gnetized region as di scussed 
bv iRempel et al.l l)2000() and. iRempel fc Dikpatil l|2003j) . 
but not to a flow pattern with two opposite flows on 
both sides of the active region belt. 

In order to address the influence of magnetic pressure, 
we show in Eq. (|13() the cj) component of the vorticity 
equation considering a stationary flow and neglecting vis- 
cous stresses (we also neglect here the magnetic tension 
force for the reason mentioned above): 



r sm 



dz 



g d 
do 



Pra 
PO 



(13) 



Note that this equation addresses the deviation of the full 
differential rotation from Taylor-Proudman state. There- 
fore si also contains entropy perturbations associated 
with the reference state differential rotation. Writing 
Q = + E^nd Si — Sr + St, where the subscript 'r' 
refers to the reference state and the subscript 't' to the 
torsional oscillation part, we get in leading order: 



2rsin6' fl 



dz 



a. 



dz 



77- 80 ' 



Em ) (14) 

PO 



Here the first term on the left hand side is the dom- 
inant one (typically at least a factor of 5 larger than 
the second one for the examples discussed here), relating 
directly the z derivative of fit to the entropy perturba- 
tion and magnetic pressure. Therefore, if the right hand 
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side of Eq. H14|l is zero, the torsional oscillation pat- 
tern has to be very close to the Taylor-Proudman state 
with d^t/dz = 0. The occurrence of St together with 
Pm/po nieans that a toroidal magnetic field with neutral 
magnetic buoyancy does not lead to any deviation from 
the Taylor-Proudman state. If magnetic buoyancy is not 
compensated, magnetic pressure causes a flow pattern 
with two opposite flows on both sides of the active region 
belt; however, they would have the wrong sign (fast mov- 
ing band on poleward side, slow moving on equatorward 
side; see also the discussion in next paragraph) when we 
make the reasonable assumption that the active region 
belt has a higher magnetic pressure than the surround- 
ing region. We will show in the following paragraph that 
a reasonable flow perturbation results from the assump- 
tion that the active region belt has lower gas pressure, 
which can be the consequence of a thermal perturbation. 

The zonal flow presented in Figure EJd) shows the 25° 
inclination angle imposed through the thermal forcing 
function. Even though the forcing function is symmetric 
with respect to the latitude of the maximum (given by 
9 — C(''))i the resulting zonal flows show a significant 
asymmetry, with the equatorward faster rotating band 
more confined than the poleward slower rotating band. 
This is a consequence of the fact that Eq. 113|l considers 
the radial and the latitudinal component of the Coriolis 
force together. The relative contribution of both com- 
ponents depends on the latitude. In high latitudes we 
have d/dz '~ d/dr, resulting in a symmetric fiow pattern 
with the slower rotating band poleward and the faster ro- 
tating band equatorward, while in low latitudes we have 
d/dz ~ —d/dO leading to a solution Vlt ^ —St, showing 
only a fast rotating band centered around the cooling re- 
gion. In mid-latitudes the solution is a combination of 
both, with having a stronger fast rotating band that is 
also more close to the center of the cooling region. Ap- 
plying the thermal forcing only to the surface, as shown 
Figure Eli), produces a zonal fiow perturbation that is 
much more confined to the surface; the 25° inclination is 
less visible in this case. 

This study shows that it is much easier to produce a 
zonal fiow pattern with the observed properties of the 
low latitude torsional oscillations through thermal forc- 
ing rather than mechanical forcing, but it does not ex- 
plain the observed 25° inclination angle, since this angle 
has to be imposed on the thermal forcing function. This 
could be a consequence of anisotropic thermal heat dif- 
fusivity that is l arger parallel than perp endicular to the 
axis of rotation l|Kitchatinov et anil994f) . A thermal sig - 
nal imposed at the surface, as sugg ested bv lSoruiti l|2003j) , 
would therefore not penetrate radially into the convec- 
tion zone but rather in an angle more aligned with the 
rotation axis. Another possible explanation could be a 
combination of mechanical and thermal forcing. To illus- 
trate this, we computed in FigureHa zonal flow resulting 
from both types of forcing. Panel a) shows the result of 
the mechanical forcing alone, panel b) the result of the 
thermal forcing. In both cases we used tq — 0.9 Rq and 
9o = 25°. The forcing functions show a radial alignment 
(C(r) — Oq — const.). The amplitude of the zonal flow 
in a) and b) is roughly the same. Applying mechani- 
cal and thermal forcing together results in a zonal flow 
with roughly 25° inclination with respect to the axis of 
rotation. 



In this section we focused entirely on static forcing 
functions, not addressing the equatorward propagation 
of the signal in the course of the solar cycle. This could 
be easily addressed by incorporating a time dependent 
phase into Eq. HlOfl : however, this does not impact the 
conclusions of this section concerning the importance of 
the Taylor-Proudman theorem, since the timescale for 
the establishment of the Taylor-Proudman state is sig- 
niflcantly shorter than the solar cycle. 

5. DISCUSSION 

We present in this paper a series of numerical exper- 
iments addressing the question of what type of forcing 
is required to produce a zonal flow variation showing 
the properties of the observed solar torsional oscillations. 
Since this study is independent from a speciflc dynamo 
model, we can infer a few general conclusions about tor- 
sional oscillations and their possible origin: 

• The Taylor-Proudman theorem applies also to per- 
turbations of ri. Therefore it is crucial to solve 
in models addressing the origin of torsional oscilla- 
tions the full momentum equation, including vari- 
ations of the meridional fiow. 

• Mechanically forced zonal flows always have the 
tendency to spread through the entire convection 
zone parallel to the axis of rotation. 

• Significant deviations from the Taylor-Proudman 
state (as observed in low latitudes) require thermal 
perturbations. A pure mechanical origin of the low 
latitude torsional oscillations is very unlikely. 

• The fact that the phase of the torsional oscillation 
signal in low latitudes is inclined such that the sig- 
nal occurs for a fixed latitude at the base of the 
convection around 2 years prior to the surface sig- 
nal does not indicate necessarily that the origin has 
to be at the base of the convection zone. This is a 
consequence of the infiuence of rotation. 

• The high latitude branch of torsional oscillations 
does not require a forcing propagating poleward in 
the course of the solar cycle. A poleward propagat- 
ing zonal fiow pattern is the natural response of the 
coupled differential rotation / meridional fiow sys- 
tem to a periodic forcing (mechanical or thermal) 
in mid-latitudes. 

• The poleward propagation speed of the torsional 
oscillation pattern is only determined by the value 
of the turbulent viscosity of the reference model. 
Best agreement is found for values of ft less than 
3 X lO^m^s^i. 

In this paper we did not consider the physical origin 
of the imposed mechanical or thermal forcings. Possible 
physical explanations for the forcings considered in this 
paper are the following: 

Macroscopic Lorentz force: 

Recently Rcnipcl (2006) showed that the macroscopic 
Lorentz force in a flux transport dynamo leads to a pe- 
riodic source in mid-latitudes, where the poloidal mag- 
netic field is sheared by the differential rotation while 
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getting transported downward by the meridional flow. 
The model provides a poleward propagating branch with 
the correct amplitude and phase relation to the magnetic 
cycle. Unlike some of the processes listed below, this type 
of feedback is unavoidable for a dynamo that is energeti- 
cally driven through differential rotation. Rempel (JOQJ) 
showed that the amplitude of the solar torsional oscilla- 
tions is consistent with a dynamo that converts around 
0.1% Lq of energy and produces around 15 kG of toroidal 
field at the base of the convection zone. 

A differe n t model w as presented recently by 
iCovas et all l|2?M l20nl [2005). Their model uses 
the observed solar rotation profile and allows for 
perturbations around this profile by considering the 
longitudinal component of the macroscopic Lorentz 
force, while the azimuthal component of the momentum 
equation is neglected (no meridional flow is considered). 
Their results show torsional oscillations with both 
branches driven by the Lorentz force that have a re- 
mark able agreement with observations f Vorontsov et al.l 
12003ft . Even though most of the magnetic activity is 
in their model concentrated in low latitudes, there are 
also magnetic field patterns moving poleward together 
with the torsional oscillation pattern. We showed in this 
paper that a poleward propagating magnetic pattern 
is not required to drive the high latitude torsional 
oscillations, but on the other hand we cannot rule out 
such a possibility based on the results we presented 
here. It would require a more detailed analysis of their 
model to determine whether the poleward propagating 
magnetic activity is the main driver or a respo nse to a 
period ic forcing in mid-latitudes as shown by iRempell 

As far as it concerns the low latitude branch of tor- 
sional oscillations we expect a significant influence of 
rotation, which is not present in a model not consid- 
ering self-consistently the meridional flow. The align- 
ment of the phase with the rotation axis makes it dif- 
ficult to obtain a surface pattern of torsional oscilla- 
tions propagating all the way to the equator especially 
if the Lorentz force action takes place in the lower 
part of the convection zone (which is the case in a dy- 
namo model living mainly on the radial shear in the 
tachocline) . This situation is comparable to the develop- 
ment of mean field models for differential rotation more 
than a decade ago. While early models that were not 
considering the meridional flow got a remarkable agree- 
ment with observations, later models solving all compo- 
nents of the momentum equation ended up in the Taylor- 
Proudman state. The problem is still under discussion, 
even though most people agree now that thermal per- 
turbations explain the obse rved profile o f the sol ar dif- 
feren t ial rotation (Kitchatinov &: RiidigeriJ995: Rempell 
120051 IMiesch et al.ll2006|) . In our experience it is much 
more difficult to explain the low latitude torsional oscilla- 
tions when the azimuthal components of the momentum 
equation are considered. The situation is different for 
the high latitude branch, since there the alignment of 
perturbations with the axis of rotation does not impose 
such a strong constraint. 



E nhance d radiative loss in active region belt: 

iSpr^iid (HOOS) proposed that the torsional oscillations 
are a response to enhanced radiative losses in the active 
region belt du e to small scale magnetic flux elements. 
iRempe f (2006) parametrized this idea and showed that 
the cooling of the active region belt drives a zonal flow 
that is, at least close to the surface, in good agreement 
with observations. Also meridional flow changes, corre- 
sponding to an inflow into the active region belt, reason- 
ably agree with observations (Komm et alJll993: iKomml 
.1994: .Zhao & Kosovichev 2004). However, the flow pat- 
tern deeper in the convection zone does not show the 
correct phase relation and depth dependence, which is 
most likely a consequence of the simple diffusive treat- 
ment of thermal perturbations. Better agreement could 
be achieved if the thermal perturbation imposed in the 
surface layers penetrates deep enough and is also influ- 
enced by anisotropic diffusivity to yield the observed in- 
clinations with respect to the axis of rotation. 

Quenching of convective energy flux: 

Due to the very small amplitude of temperature per- 
turbations required for driving zonal flows it seems con- 
ceivable that quenching of convective energy flux by the 
dynamo generated field can also cause thermal shadows 
within the convection zone resultin g in zonal fiow s. Con- 
trary to the process proposed by ISpruil l)2003j) . these 
temperature perturbations would originate close to the 
base of the convection zone, where the strongest mag- 
netic field is found. Again, an anisotropic convective 
energy flux would be required to explain the observed 
phase relation. 

Quenching of turbulent angular momentum flux (A- 
quenching): 

Microscopic Lorentz force feedback in terms of quench- 
ing of turbulent transport proces ses driving differen- 
tial ro tat ion has been addre ss ed by iKitchatinov fc Pipiii 
(119981), I Kitchatinov et alJ (|1999D ' and iKiiker et al. 
(1999) as possible explanation for torsional oscillations. 
Since A-quenching leads to an additional source term in 
the equation for f2, this is a mechanical forcing similar 
to macroscopic Lorentz force feedback. Therefore all the 
results discussed above for mechanical forcing also apply 
in this case. 

It has been speculated whether the similar inclination 
of the phase of the torsional oscillation pattern and the 
isorotation contours is a coincidence or not. None of 
the processes discussed here would necessarily lead to a 
similar angle; however, the physical cause behind both, a 
thermally induced deviation from the Taylor-Proudman 
state, is similar. 



This work originated from discussions with Manfred 
Schiissler during a visit of the author at the Max-Planck 
Institute for Solar-System-Research in Katlenburg- 
Lindau, Germany. Very helpful comments by M. 
Schiissler, K. B. MacGregor, P. A. Gilman and the 
anonymous referee are gratefully acknowledged. 
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